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The two-body density matrix for "^He,^^ O and ^^Ca within the Low-order approximation of the 
Jastrow correlation method is considered. Closed analytical expressions for the two-body density 
matrix, the center of mass and relative local densities and momentum distributions are presented. 
The effects of the short-range correlations on the two-body nuclear characteristics are investigated. 
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Nowadays two-body knock-out reactions such as (j,NN) [|l|j2| and {e,e'NN) are intensively studied in order 
to extract some information about the nucleon-nucleon short-range correlations (SRC) in nuclei. In particular, the 
cross section for the two-nucleon emission processes is generally related to the two-body spectral function at least 
in PWIA 0. The first calculations of the two-nucleon spectral function of ^^O have already been performed Q by 
O i' treating the long-range correlations within a Dressed RPA approach and including the SRC in terms of the defect 
functions emerging as a solution of the Bethe-Goldstone equation for finite nuclei. The results have been successfully 
applied for calculating the (e, e'pp) cross-sections Q. 

Due to the complexity of the problem, however, it is highly desirable to reach a theoretical description of the 
' two-body knock-out reactions directly in terms of the nuclear ground state (e.g., the ground-state two-body density 
matrix) without the necessity to deal with the two-body spectral function which is an enormously more complicated 
object due to the presence of various excited states of the system. Making use of series of more or less controlled 
approximations, people usually try to incorporate in this context simplified expressions or combinations of physical 
t — I quantities such as the two-particle spectroscopic factors and overlap functions j^, relative and center-of-mass pair 
' momentum distributions, combined two-body momentum distributions pc| ] and generalized momentum distributions 
"nI (see e.g. [0). Realistic but still elaborated description of the ground state properties of nuclei up to ^^O is known 
from the Variational Monte-Carlo calculations and recently from the Green's Function Monte-Carlo method 

performed for nuclei with mass number A < 7 [|15|| . In any case, one still needs a simple and possibly analytical model 
description of the ground state two-body density matrix in order to clarify the two-body nuclear characteristics as 
^ I well as the extent to which they are affected by the nucleon-nucleon correlations. 
I . The situation is quite similar to the theoretical description of one-particle removal processes in terms of the nuclear 
Q ' ground state characteristics. At the beginning, it has been demonstrated [|6| that the knowledge of the ground-state 
^ , one-body density matrix of the target nucleus is sufficient to restore the single-particle overlap functions, spectroscopic 
' factors and separation energies associated with the bound {A — l)-particle eigenstates. Then, quantitative estimates 
^ [ have been obtained ||l7[| within a simple and analytical one-body density matrix model [l^ which takes into account 
the short-range correlations in terms of the Jastrow correlation method. Plausible conclusions have been made 
for the properties of singe-particle overlap functions in comparison with the associated shell model orbitals and 



^ I the natural orbitals [^^^ which are of frequent interest in this context |^,^. The resulting overlap functions 
■ " ' and spectroscopic factors have been used to analyze differential cross sections of {p, d) reactions and single-particle 
momentum distributions in {e,e'p) reactions p3[ . Finally, more sophisticated representations of the one-body density 
matrix [p^-p6| have been used for extracting overlap functions and spectroscopic factors further applied for analyzing 
{p,d) reaction cross-sections |2^]. Thus, the resulting comparative study j27j has clarified the impact of the different 
types of nucleon-nucleon correlations on the one-particle removal reaction cross-sections. 

Recently, similar restoration procedure has been proposed in Q connecting the two-nucleon overlap functions 
associated with the bound states of the {A — 2) particle system with the asymptotic behavior of the ground state 
two-body density matrix of the A-particle system. This makes it possible to calculate, at least in principle, the two- 
body overlap functions, spectroscopic factors and separation energies on the basis of realistic representations for the 
ground state two-body density matrix. Again, as a first step, one clearly needs a simple and analytical representation 
of the two-body density matrix which adequately reflects the properties of the realistic nuclear ground state and also 
takes into account the short-range correlations in nuclei. 

The present paper suggests such a simple and analytical representation of the ground state two-body density matrix 
derived within the Low-order approximation of the Jastrow correlation method. The associated two-particle nuclear 
characteristics for the closed s — d shell nuclei '^He^^^O and ■*°Ca are analyzed. A comparison with the realistic 
Variational Monte-Carlo calculations is also made. It justifies the physical meaning of the approximation as a tool 
for analyzing the impact of the nuclear SRC in a simple way. The analytical expressions obtained can facilitate the 
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explanation of most of the questions arising from the interpretation of physical properties which are important in 
treating the two-particle emission processes in nuclei. 

The paper is organized as follows. The general definitions of the two-body density matrix and the associated nu- 
clear characteristics in coordinate and momentum representation arc introduced in Sec. II. The analytical expressions 
derived within the Low-order approximation of the Jastrow correlation method are given in Sec. III. Closed analytical 
expressions for the two-body nucleon momentum and density distributions are collected in Sec. IV. Results and dis- 
cussions are given in Sec.V, while the conclusions are summarized in Sec. VI. The Appendix contains the coefficients 
entering the analytical expressions for ^TJe and ^^O nuclei. 



II. TWO-BODY DENSITY MATRIX 

A. Definitions and properties 

The physical antisymmetric state of a system of A identical fermions 'i'^ normalized to unity defines a set of density 
matrices of order p = 1,2, A 

p^P\xi,X2, ...,Xp;x[,x'2, ■.■,x'p) =< '^{a^ {xi)a\x2):-a^ {xp)a{x[)a{x2)...a{x'p)\'^ > , (1) 

where a^{xi) and 0(2;^) stand for creation and annihilation operators for a nucleon at position Xi, which includes the 
spatial coordinate r^, the spin Si and the isospin r^. In particular, the one- and two-body density matrices are defined 
in coordinate space as: 

p'^^\x,x') = {■9'^^^\a+{x)a{x')\^^^^) , (2) 

and 

p^^\x,,X2;x[,x',) = {¥''^\a+{xi)a+ix2)aix[)a{x',)\¥^^) , (3) 

respectively. From these defining equations one can easily recognize many of the properties of the density matrices. 
They are Hermitian 

p^^\x,x') = p'^^^* {x ; x) , p'^^Hxi,X2;x[,x'2) = p'^^^*{x[,x'2;xi,X2) , (4) 
and trace-normalized to the number of particles and of pairs of particles: 

TVpti' = / p{x)dx ^ A , (5) 



Trp(^) = 1 1 p^^Hx,,X2)dx,dx2 = ^^^^ , (6) 
with diagonal symmetric elements 

p{x) ^ p''^^x,x) , p'^^\xi,X2) = p^'^\xiX2;XiX2) . (7) 

In addition, the two-body density matrix p^^' is antisymmetric in each set of indices, e.g., 

p'^^\xiX2;x[x2) = -p^^\x2Xi;x[x2) , (8) 

so that its diagonal elements vanish identically if both coordinates are equal, i.e., p^^'^ {xi, Xi) = 0. 
The one- and two-body density matrices are related by the formula 

p^'^\xiX2; x[x2) dx2 = ~~Y~ '"^^''(^I'^i) ' (9) 
and both can be presented in the momentum space using the Fourier transforms: 

n^^\k;k') = j p{x,x') exp (k.r - k'.r')] drdr' (10) 



(fci, fc2; /c'l, /C2) = y p'-^-'(a;i,2:2;xi, a;2) exp [i (ki.ri -I- k2.r2 - k'l.r'i - kj.rj)] dridr2dr\dr'2 , (11) 

where ki stands for the momentum k^, spin Si and isospin of the z— th particle. Relations similar to eqs.(^)-(^) held 
for the one- and two-body density matrices in the momentum space as well. 
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B. Two-body nuclear characteristics 



Typical ground state quantities of interest one usually considers arc the local density 

p{r) = p{r,r) = ^p{x,x), (12) 

ST 

the associated elastic form-factor 

i^(q) = jj p{r) exp[iq.r] dr (13) 

and the nucleon momentum distribution 

n(k) = n(k, k) = ^ n{k,k) (14) 

ST 

obtained after spin (and/or isospin) summation of the diagonal elements of the one-body density matrix in coordinate 
and momentum representation. 

The two-particle emission experiments however require some knowledge of physical quantities associated with the 
two-body density matrix. For example, the diagonal elements of the two-body density matrix p^^' in coordinate space, 
eq.(|^), define the center-of-mass pair local density distribution: 

p(2)(R) = J /2)(R + s/2,R-s/2)dr (15) 

and the relative local density distribution: 

P^^\s)= j p(2)(R + s/2,R-s/2)rfR (16) 

while the diagonal elements in momentum space eg. (pi]) define the associated center-of-mass and relative pair mo- 
mentum distributions: 

n(2)(K)= j ?i(2)(K/2+k,K/2-k)dk (17) 

and 

„(2) (k) = j (K/2+k, K/2-k)dK , (18) 

respectively. 

The physical meaning of p^^^ (s) and rt^^^ (k) is the probability to find two particles displaced of a certain relative 
distance s = ri— r2 or moving with relative momentum k = (k-^— k2)/2, respectively, while p*^2^(R) and rS'^'^^K.) rep- 
resents the probability to find a pair of particles with center-of-mass coordinate R = (ri + r2) /2 or center-of-mass 
momentum K = ki-t-k2, respectively. 

III. ANALYTICAL EXPRESSIONS FOR THE TWO-BODY DENSITY MATRIX 

A. Mean-field approximation 

The mean-field approximation to the nuclear ground state of an A particle system is represented by a single Slater 
determinant 

^sd{xi,X2,---,xa) = -^dei\'fi{xj)\ , (19) 

where the orthonormalized set of single-particle functions ipi{x) = ipi{r,s,T) is emerging from some kind of shell 
model or self-consistent mean-field calculations. The ground state ^sd incorporates two kinds of correlations: (1) 
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Pauli correlations associated with the antisymmetric properties of ^gjj and (2) the correlations among the nucleons 
forming the nuclear mean field that determines the particular form of the single particle states fi{x). 

The following expressions for the one- and two-body density matrices are well known from the mean-field theory 

A 

Psd(x,x') , (20) 



A 

.(2) 



« J = l 

where the antisymmetric uncorrelated two-body wave functions are used 

ipij{xiX2) = [ip*ixi)'P3ix2) - (Pj{xi)(prix2)] ■ (22) 

(2) 

In order to obtain analytical expressions which will allow us to compute psn and pgjj in a direct way we further 
assume all states belonging to the uncorrelated Fermi sea as represented by harmonic oscillator single-particle wave 
functions Lpf^{x) which depend on the harmonic-oscillator length a, heaving the same values for protons and neutrons. 
In particular, these are the states Is for "^i/e, Is and Ip for ^^O and Is, Ip, Id and 2s for '^^Ca. Because we are interested 
in spin (and/or isospin) free quantities like total center of mass and relative coordinate and momentum distributions 
we consider only the matrix elements which are fully diagonal in spin and isospin variables. Under these assumptions 
eq.(|2|) has closed analytical form: 

PsD(i"i''5i''^i;''2,S2,r2;r3,si,Ti;r4,S2,r2) = - [p(ri, r3)p(r2, r4) - (5rir2^siS2P(ri, r2)p(r3, r4)], (23) 
since the spin and isospin free one-body density matrix 

p(ri,r2)= ^ps_D(i-i,s,r;r2,s,r) (24) 

ST 

is an explicit product of exponent and polynomial factors depending only on the scaled coordinates — ar^: 

Q/3 ~|~ X^ 

P^^i-'^^) = exp[ ^ ^ ^ ] PsD (xi,X2) , (25) 

where 

PsD (xi,X2) = 1 for "ffe , 

= (1 + 2x12) for , (26) 

= i [5 + 42)12-2 (a;? -2x12+ a;2)] for ^^Ca 
and Xij = x^.Xj = a^rirj cos 9ij, 9ij being the angle between radius vectors r,j and r, . 



B. Jastrow correlations method 

In the present paper we consider the nucleon-nucleon SRC within the Jastrow correlation method p9|-p2|. This 
method incorporates the nucleon-nucleon SRC in terms of the wave function ansatz: 

¥^Hrur2,...,rA) = iCA)-'/' J] /(I r. - r, |) $^^(ri, r2, . . . , r^), (27) 

l<i<j<A 

where ^gjj is a single Slater determinant, /(r) is a correlation factor which goes to unity for large values of r and Ca is 
a normalization constant. Except for systems containing very small number of particles it is impossible to calculate the 
one- and two-body density matrices using the Jastrow ansatz eq.(p7|). One usually applies a perturbation expansion 
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in terms of linked diagrams |Q . The so-called Low-order approximation (LOA) keeps all terms up to second order 
in /i = / — 1 and first order in g = — 1 in such a way that the normalization of the density matrices is ensured 



(28) 



order by order |gO|. In particular, the resulting two-body density matrix is of the form pi 32|: 

p^PoAi^,,X2;x3, X,) = p^'i,(1234) + [ f^, - 1] ^^'^,(1234) + / dx^ [{f^.^^ - 1) + {1^45 - 1)] 

('21 ('21 (2.) 

PSD {xi ; a;3)p^f, (2545) + pso {xi ; X4)p^g^ (2553) -I- pso {xi ; X5)p^g^ (2534) 
+jjdx^dx^ [/5*6/56 - 1] {pi'i,(2665)p^'i,(1534) + PsD{xi;x5)psD{x2;x:!,)p''gl,{bMQ) 

('2) ('21 1 

+ PSD (xi ; x^)psD {x2 ;x4)psD (5663) + PSD {xi ; X5)psD {x2 ; xe)p^gly (5634) j , 

where = /(|rj — fjl), Pg£)(1234) = ^^^(xi, 0:2; xs, 0:4) and the integration over Xi means summation over the spin 
and isospin variables and integration over the spacial coordinate. 

It is clearly seen from eq.(|2^) that P^j^qj^ generally depends only on two ingredients, the correlation function 
/(|ri — Tjl) and the two-body density matrix in its mean- filed approximation , eq. (|2l]). In the previous 
eqs.( |25|p^ ) we have already derived in closed analytical form in terms of harmonic-oscillator single-particle wave 
functions. In order to find closed analytical expression for p^^q^ , eq.(|2|), we further assume that the correlation 
factor /(r) is state- independent and has simple gaussian form: 

/(r) = fir) =1-0 exp(-/3V2) , (29) 

where the correlation parameter P controls the healing distance, while the parameter c accounts for the strength of 

the SRC. Under these additional assumptions, performing explicitly the integrations entering eq.(p8|), the diagonal 

(2) 

part in the spin and isospin variables of p\qj^ transforms to pure algebraic expression: 

(2) / X (2) , (2) , (2) , (2) ,on\ 

where the first term p|j^^ is already defined by eq.( p5| , p6| ). The sum of this term p^^^ with the next one 
Pa = 1 {c^ea;p[-(zi - Z2)^ - (Z3 - Z4)2] - cea;p[-(zi - Z2)^] - cea;p[-(z3 - •Lif]) 

(31) 

X {p(ri,r3)p(r2,r4) - (5^i^2(5siS2p(ri,r2)p(r3,r4)} , 

where = I3vi is often referred to as a first order approximation to the Jastrow two-body density matrix ]3^ . The 
main disadvantage of this approximation is that it does not satisfy the normalization condition cq.(|^) and thus has 
restricted physical significance. The third term p^"* in eq.(|30|) has the form: 

P^B = -r {4 [p(ri, ra) p{^2, ^i) - (5rir2'5siS2/o(ri, i"4) p(r2, ra)] (-^13 + hi) - p{^i,^3.) (-^1324 + ^2424) 

- p(r2, r4) (/i3i3 + I2) + (5rir2'5siS2 [p(r2, r3) (-^1314 + ^2414) + P(ri, I"3) (-^1324 + -^2424)]} 

-f {4[p(ri,r3)p(r2,r4) - 5r^r25siS2 p(ri, r4) /o(r2, r3)] {h +h +h + h) (32) 

-P(ri, r3) (/i24 + ^224 + ^324 + -^424) - p(r2, r4) (/ll3 + ^213 + ^331 + -^413) 

[p(r2, r3) (/ii4 + /214 + ^314 + ^414) - p(ri, r4) (/123 + /223 + ^323 + -^423)]} , 

(2'] 

while the last term p)j' reads: 

p[?) = cA^(ri,si,ri;r2,S2,T2;r3,si,ri;r4,S2,T2; y) ~ \ c2A^(ri, si, ti; r2, S2, T2; r3, si, n; r4, S2, T2; 2y) , (33) 



M =4[p(ri,r3) F24-p(r2,r4) Y13] - Y1324 + p(ri, r3) ^42 -p(r2,r4) Z13 

-<5riT2'5siS2 {4 [p(ri,r4)F23 - p(r2,r3)Yi4] - 11423 - p(r2,r3)Zi4 - p(ri, r4)2'23} • 



(34) 
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where the following products of exponential and polynomial factors depending on the scaled coordinates = aVi and 
the dimensionless parameter y = j o? are introduced: 



and 



h 
hi 

-fl234 



eM-Tl^]Pii (xi) 



exp[- 

^3/2 ^^P 



^ 1+2^ ^JP/2(X1,X2) 



(35) 



+ + y (2a; ^+3;^ +3:3) 1 p ( \ 
2{l+y) i -f/3lXljX2,X3j 



exp[ 



X3^ + Xj'^ + 2y {xi^+X2^) + 2y''(xi-X2)^ 
2(l+2ij) 



] f/4(xi,X2,X3,X4) 



>"i2 = exp[ 
Z12 = ^ exp[ 
^1234 = ea:p[- 



]Py2(xi,X2; y) 



-^^^4^]P22(xi,x2; y) 



(36) 



Fy3(xi,X2,X3,X4; y) 



The two-body density matrix represented by eqs.(|30|-p6D does not depend on the particular choice of the nucleus. 
All nuclear structure information is reflected by the polynomials {Pi, Py, Pz} entering eqs.(^ 36). Their explicit form 
is given in the Appendix for ^iJe and ^^O nuclei. As a result we end up with closed analytical algebraic expressions 
for the two-body density matrix involving three independent parameters, the harmonic oscillator length a and the 
correlations parameters c and /3. Due to its gaussian-polynomial structure one can easily derive the Fourier transforms 
defining the associated two-body density matrix in the momentum space. 



IV. LOCAL DENSITY AND MOMENTUM PAIR DISTRIBUTIONS 

Having derived the two-body density matrix in the coordinate space one is able to calculate all two-body nuclear 
characteristics in close analytical form. Changing the coordinates of two particles ri and r2 in eqs.( ^0[]3^ ) to the 
center-of-mass R = (ri + /2 and relative s = (ri— r2) coordinates, one obtains the center-of-mass ^^^•'(R) , eq.(^5|), 
and relative p^'^\s) , eg. jl^) , pair local distributions. Using the dimensionless variables R = ajRI and s — a\s\ they 
have the following form: 



p'-'^^R) ^ Aexp -{V2rY +Bexp 



2{l+2y) 
2+3 y 



2R) 



-C exp 



1+4 y 

1+3 y 



2R) 



(37) 



where 



p^^\s)=Vexp -{s/V2) +£exp 



2(1+2;/) 
2+3 y 



'-Texp 



-f^tf(^/V2)^ 



+Qexp 



-(l + 4y) (s/V2)' 



— Ti, exp 



-{l + 2y) {s/^y 



A^^['m{R) + 2cv2{R,y)~ rj2{R,2y)], B ^ ^crj^iR), C= -^c^f^^{R), 



(38) 



(39) 



= ^ [mi(s) + 2cfi2is;y) - fi2is;2y)] , £ = ~:^2cp3{s;y), T = -^c^ /i3(s;2y), 

a-^cVi(s), W = -^2c/ii(s). 



(40) 



The center-of-mass n(^)(K) , eq. (p^) , and relative ?i'^^^(k) , eq.(p^), pair momentum distributions follow after 
performing Fourier transform of the center-of-mass p(^) (R,R') and relative p^^) (g^g') pair density matrices, respectively. 
Again, introducing K = |K| ja and fc = | k| ja one obtains: 
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n(2)(K) ^Mexp 



B' exp 



1+2 y 
2+5 



y (^^) 



-C exp 



(41) 



where 



2(1+2 y) 
2+5 y 



(k) = P' eip - (\/2fc) + £' exp 

+G'exp {V2kf] -H'exp -i±|f (V2fc) 



-4' = 5^ (71 (^) + 2c72(X, y) - c2 72 (if, 2 y)) , = ^-^c73(i^) , C'= ^^c^ 74 (if), 



(42) 



(43) 



V 



2q3 



475(0i(fc) + 2c02(fc,2;)- c2 02(fc,2y)), £'= c04(fc), ^' = ^^^c^ (fc), 



Q' - 211^^2 05 (fc), W = ^^ce,{k) . 



(44) 



In the above expressions the exponential dependence is expUcit while the expressions for the polynomial amplitudes 
{r]i}, {ni}, {ji} and {Ot} are presented in the Appendix for '^He and ^^O nuclei. 

For completeness, we are giving also the local nuclear momentum distribution n{k), eq.(p^), which is associated 
with the one-body density matrix. Its closed analytical form within the present model has already been derived in 



i(k) = A" exp + S"exp 



1 



l + 2t/ 



- C'cxp 



1 + 22/ .2 
1 + 3 y 



(45) 



and the polynomial expressions for the amplitudes A",B",C" can be found in Q. Considering the gaussian factors 
one can realize a factor V2 which scales the momenta entering the three momentum distributions n'^^)(K), n(k) and 
n(2)(k). 



V. RESULTS AND DISCUSSION 



Next step in our study is to define the parameters of the problem, the oscillator parameter a and the parameters /3 
and c related to the healing distance and the strength of the SRC, respectively. The rigorous procedure would be to 
apply the variational approach based on the one- and two-body density matrices defined so far by minimizing the total 
energy of the system with a realistic A^A'^— interaction with respect to the parameters a, P and c. As in our previous 
papers p^], however, we prefer to obtain the values of a and P phenomenologically by fitting the experimental elastic 
formfactor data using the analytical expression for the elastic formfactor following within the present model from 
eq.(|l^). The values of the parameter c for '^He and ^^O are determined under the additional condition the relative 
pair density distribution p^^^ (s) , eq . (^8|) , to reproduce at s = the associated value obtained within the Variational 
Monte-Carlo approach [|l3[ and [ |l4| , respectively, while we simply put c = 1 for ^^Ca since there are no variational 
calculations for '^^Ca . Thus, in the present numerical calculations we are using the following values of the parameters: 



a = 0.82 fm 
a — 0.61 fm 
a = 0.52 fm 



-1 

1 
1 



(3 = 1.23 fm-\ 
P = 1.30 fm-\ 
P = 1.21 fm-\ 



c = 0.76 for ^iie , 
c = 0.77 for , 
c = 1.00 for 40(70 . 



(46) 



At the beginning we are comparing in Fig. 1 some results from our crude but analytical model (the solid lines) 
with the results (dashed lines) emerging from the orders of magnitude more complicated Variational Monte-Carlo 
calculations ]l3[|, [T^ . The top panel of Fig. 1 clearly demonstrates that our simple gaussian form for the correlation 
factor /(r) , eq. (|29[) , should be considered only as a first approximation to the actual shape of the realistic central 
correlation function. Nevertheless, the analytical results for the relative density pair distributions in '^He and ^^O ( 
two bottom panels in Fig. 1) show quite acceptable quantitative agreement with the Variational Monte-Carlo results. 
Despite of the simplicity of the model and the LOA used, the characteristic behavior of the realistic pair density is 
obviously reproduced. The small values of the relative pair distributions at s = indicate the presence of significant 
SRC which in the Variational Monte-Carlo calculations are due to the repulsive core of the two-body interaction. 
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In Fig. 2 a comparison is made between the correlated (c 7^ , solid curves) and the uncorrelated (c = 0, 
dashed curves) results for the relative pair density distribution p*-^'(s) , eg. (|3^) , for *He, ^^O and ^"Ca. It is 
seen that due to the SRC the shape of the distributions changes significantly. The SRC lead to a deep hole in 
the correlated distributions near s = 0, while the uncorrelated distributions saturate at small distances to values 
which are significantly different from zero. Our calculations have shown that the SRC do not affect significantly the 
center-of-mass pair local distributions in the coordinate space. 

The above examples show that despite of the simplicity of the model it is able to incorporate the SRC into the 
two-body quantities of interest. Therefore it may be considered as a useful starting point in analyzing two-particle 
emission experiments. 

From (e, e'2N) experiments one can expect important information about the effects of SRC on the pair local 
distributions in the momentum space. In Fig. 3 we present the center-of-mass (solid lines) n^'^\K), eq.(pl|), and 
relative (dot-dashed lines) ri(^^(fc), eq.(^), pair momentum distributions for '^He, ^^O and '^'^Ca . A comparison is 
also made between the correlated (c ^ 0) and uncorrelated (c = 0) results. For completeness, by dashed lines in Fig. 
3, we also show the results for the local nuclear momentum distribution n(fc), eq.(|5|), which is associated with the 
one-body density matrix. 

From Fig. 3 is clearly seen that due to the SRC high momentum tails develop at large values of the momenta and this 
behavior is typical for all three kinds of momentum distributions considered. Comparing with the uncorrelated results 
one can see that the correlation effects start to dominate first for the relative pair momentum distribution rS'^\k). 
Then, at larger values of the momenta, the center-of-mass distribution rS^\K) also develops a high momentum tail. 
The local nuclear momentum distribution nik) takes an intermediate position between both distributions rS^^k) and 
rS'^\K). In the case of '^He for example the SRC become important for momenta larger than 1.4, 2.1, 2.7 fm~^, 
for the momentum distributions n'^'(fc), n(fc) and n(^)(i^), respectively. Similar values are valid also for nuclei ^^O 
and ^°Ca. The reason for such behavior is obviously due to the scaling factor of 1/2 which has been observed in 
the analytical expressions for n'^'^\K), eq.(|4l|), n{k), eq.(p]) and n^'^\k), eq. (^^ . Also, at large momenta both the 
correlated and uncorrelated momentum distributions obviously satisfy the inequality: 

n(2)(fc) < n{k) < n^'^\K) . 

Of course, the comparison of the results given in Fig. 3 has to be done keeping in mind the different meaning of 
the arguments of the three types of momentum distributions considered.. 

Plotting together each of the nucleon momentum distributions n'^^^ (fc), n(fc) and n^^^^ [K) for all nuclei ^He^ ^^O and 
^^Ca, as it is done in Fig. 4, one can detect the interesting fact that in the high momentum region these distributions 
are almost universal in the sense that they do not significantly depend on the mass number A. This universal behavior 
of the local nucleon momentum distribution nik) has been observed earlier (see e.g. |p^ , p^ , p2| ) . Obviously similar 
tendency exists also for the relative and center-of-mass pair momentum distributions. 

VI. CONCLUSIONS 

In this paper we have derived closed analytical expressions for the two-body density matrix and associated two-body 
nuclear characteristics within the LOA to the Jastrow correlation method for the closed s — d shell nuclei '^He,^^0 
and '^^Ca under two simplifying assumptions: harmonic-oscillator single-particle wave functions entering the Slater 
determinant and a state-independent gaussian-like correlation function. The comparison with more realistic results 
emerging from Variational Monte-Carlo calculations has shown the usefulness of the expressions derived. They can 
be applied as a starting tool in analyzing two-particle emission experiments where important information about the 
effects of SRC on the two-body nuclear characteristics is expected. In particular, it has been shown that the SRC 
effects start to dominate in the high- momentum region of the relative n^^)(fc), local n{k) and center-of-mass n'-^'>{K) 
momentum distributions at different momentum values which are proportional to a factor of \/2 and are almost 
independent on the nucleus considered. An universal asymptotic behavior of the relative and center-of-mass pair 
momentum distributions (normalized to unity) is indicated similar to the well known high-momentum tail of the 
nucleon momentum distribution n{k). 

As a first approximation, the simplicity of the results can help us to concentrate our attention towards the com- 
plicated questions arising from the involved physical interpretation and the mechanism of the two-particle emission 
processes in nuclei. 
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VII. APPENDIX 
A. Polynomials entering the two-body density matrix 

The expressions for {P} in eqs.(^ 36) which determine the two-body density matrix are: 



1. Nucleus '^He 
-3/2 



Pi, = Pi, = (1 + 2yy 
Pi,=Pi, = {l + y)-'/' 

Py, (y) = Pz, (y) = Py, (y) = (i + 2yy'/' 







2. Nucleus ^'^O 




PiA^) 


- (l+y)V2 [4 + 


5y + y^ (1 + 2x2)] 




f/2(xi,X2) 


_ 2 To 
" (l+2y)^/^L^ 


+ 5y + y^ + +X2^ + 2x12)] 




-P/3(X1,X2,X3) 


- (l+yy/2 [1 + 


2x23 + 2y (1 + ^12 + + 2:23 ) - 


fy2 (1 + 2x12 +2x13 +4xi2 x 


P/4(xi,X2,X3,X4) 




f- 2 0:34 -1- 2 y ( 2 -1- a;i3 + xu + X23 ^ 


- X2i +2 X34 ) 



+ 47/2 ( + ^13 + ^ _|_ 2:23 + Xi3 X24 X24)] 

Py, (xi,X2; y) = (1+2^ [4 + 8x12+2/ (13 + 18x12) +2/2 (10+14x12)] 
Pz2(xi,X2; 2/) = (TTn^fT? [1 + 2x12+2/ (7 + 6xi2)+2/2 (lO+Uxia)] 
Py,(xi,X2,X3,X4; 2/) = (i+2y)V2 [1 + 2x12 + 2x34 + 4x12X34+ 

2 2/ (2 + 3 X12 + xi3 + xi4 + X23 + X24 + 3 X34 + 4 X12 X34) 

+42/2 (1 + X12 + xi3 + xi4 + X23 + 2:24 + a;34 + X14 X23 + 2:13 X24 + a;i2 X34) ] 
3. Nucleus *°Ca 

The expressions for *'^Ca have similar structure but are up to 4-th order with respect to x.y, thus being too fong 
to be presented here. They can be obtained in the form of user friendly files for Mathematica 3.0 upon request. 

B. Polynomials entering local density pair distributions 

The expressions for {77^}, {fJ-i}, {7i} and {Oi}, which determine the two-body density and momentum distributions 
eqs.(|37|-^ are: 
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1. Nucleus *He 

m{R) = 12^/2 r,2{R; y) = -j^^^ Vs{R) = j^^i^ m{R) = - j^M)^ 

fii{s) = ^ H2{s; y) = ^(1^^,' ^)3/. M3(s; y) = (2+3^)3/^ 

01 (fc) = 4 02(fc; y) = 03(fc) = (1+4^ 

04(fc) = - (2+5 "^3/2 6'5(fc) = (1+4 J/)i* ^6(fc) = (1+ y)3/2^(l+4 j^)3/2 

2. Nucleus ^^O 
r]i{R) =2^2 (I17+104i?2 + 48ij4^ 

??2(ii; 2/) = (i+2^)V2 [48i?4(28 + 60y + 35y2) 

+ 8 i?2 (364 + 848 y + 555 y^) + 3 (1092 + 2780 y + 1837 y^)] 

Vs{R) = - (2+^15/2 [l536 iiS' (1 + 3 + 2 2/2)2 

+ 192 i?4 (28 + 224 y + 749 y^ + 1352 + 1419 + 843 + 225 y^) 

+ 8 ii^ (2 + 3 yf (364 + 1762 y + 3489 + 3105 y^ + 990 

+ 9(2 + 3y)^(182 + 563y + 588y2 + 212y3) 

^4(i?) = (j^|^ [384iiS2(i + 6y + 8y2)2 

+ 48 i?^ (7 + 112 y + 749 y^ + 2704 y^ + 5676 y^ + 6744 y^ + 3600 y^^) 
+ 8 i?2 (1 + 3 yf (91 + 881 y + 3489 y^ + 6210 y^ + 3960 y'') 
+ 9(l + 3yf (91 + 563y+1176y2 + 848y3) 

'"i(^) =1^ (93 + 34s2 + 3s4) 
-"2(5; = 4 v^(i+2^)V2 [2697+ 6540 y + 4047 y2 

+ 2 (493 + 1214 y + 831 y^) + 3 (29 + 64 y + 39 y^)] 



/^3(s; y) 
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(2+3j/)i5/2 



3 (2 + 3 y)^ (434 + 1245 y + 1158 y^ + 348 y^) 

+ 2 s2 (2 + 3 yf (476 + 2370 y + 4713 y^ + 4215 y^ + 1386 y^) 
+ 12 (28 + 232 y + 813 y^ + 1554 y3 + 1733 y^ + 1083 + 297 y^') 
+ 24s62;2(i + 3y + 2y2)' 
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^) = 480(1+2,)^/- (28 + 60y + 35y2) 

+ 2K^ (364 + 872 y + 603 y^) + 3 (l092 + 2732 y + 1741 y^)] 



^3(-?^) = i5(2+5w^ [24i^S'(l + 2yr(l + 42/ + 3y2) 

- 12 (28 + 344 y + 1705 y^ + 4340 y^ + 5933 y"* + 4087 y^ + 1105 y^) 
-2K'^{2 + 5 yf (364 + 2870 y + 8263 + 10185 y^ + 4680 y^) 

- 3 (2 + 5 y)^ (546 + 3251 y + 6160 y^ + 3770 y^) 

= 120(1+.)"/- (1+4 [3 i^M7 + 46 y + 72 y^) 

+ 2 (91 + 707 y + 2250 y^ + 3238 y^ + 1604 y^) 

+ 3 (1 + y)^ (273 + 1886 y + 5060 y^ + 6424 y^ + 3472 y^) 

6ii(fc) = ^ (117+104fc2+48A;4) 

^) = 60 (1+2 .f/- [48fc4(29 + 64y + 39y2) 

+ 8 fc2 (493 + 1202 y + 807 y2) + 3 (899 + 2204 y + 1397 y^)] 
^3^= - 30(1+4^)7/2 [48A;4 + 8fc2(17 + 80y + 24y2)+3(31 + 224y + 400y2)] 

= I^(iffS^ [l536fcS2(i + 2yf (l+4y + 3y2) 

- 3 (2 + 5 y)^ (434 + 2571 y + 4800 y^ + 2870 y^) 

- 8 (2 + 5 yf 476 + 3774 y + 11043 y^ + 14075 y^ + 6780 y^) 

- 192 A;'' (28 + 344 y + 1697 y^ + 4268 y^ + 5691 y'^ + 3727 y^ + 905 y^)] 

^5(fc) = 60ITOF [48fc4 + 8fc2(17+160y + 368y2) + 3(l + 4y)'(31 + 200y + 496y2) 

^6(fc) = i5(i+,)ii/Mi+4.)^/- [48 (7 + 46 y + 72 y^) 

+ 8 fc2 (119 + 935 y + 2862 y^ + 3890 y^ + 1844 y^) 

+ 3(1 + yf (217 + 1486 y + 4236 y^ + 5944 y^ + 3472 y^)] 

5. Nucleus '^°Ca 



Again, the expressions for "^°Ca have similar structure but are too long to be presented here. They can be obtained 
in the form of user friendly files for Mathematica 3.0 upon request. 
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Figure Captions 



Fig. 1. Comparison of the present results (solid lines) with the Variational Monte-Carlo results for 
Argonne W14 potential |Q as obtained in |14|, (dashed lines). Top panel - comparison of the central 
correlation function. Next two panels - comparison of the relative pair density distributions normalized as 
47r J p^^^s)s'^ds = 1 for '^He and ^^O, respectively. 

Fig. 2. Comparison between correlated (c ^ , solid lines) and uncorrelated (c = 0, dashed lines) 
results for the relative pair density distribution (BSh normalized as 47r J p'-^'> (s) s-^ ds = 1 for "^He, ^^O and 

Fig. 3. Correlated (c ^ 0) and uncorrelated (c = 0) results for the center-of-mass n^^^ (K) (solid curve) 
and relative n^"^^ (fc) (dashed-dot curve) pair momentum distributions and for the nucleon momentum 
distribution n(fc) (dashed curve) for '^He (top panel), ^^O (middle panel) and ^'^Ca (bottom panel). All 
distributions are normalized to unity as e.g. 47r / n{k)k'^dk = 1. 

Fig. 4. Correlated momentum distributions from Fig. 3, n{k) (top panel), n(^^(if) (middle panel) 
and n^'^\k) (bottom panel), but collected for all three nuclei ^He (solid curve), ^^O (dashed curve) and 
^°Ca (dashed-dot curve). 
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